Excitations of the bimodal Ising spin glass on the brickwork lattice 
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An exact algorithm is used to investigate the distributions of the degeneracies of low-energy excited 
states for the bimodal Ising spin glass on the brickwork lattice. Since the distributions are extreme 
and do not self-average, we base our conclusions on the most likely values of the degeneracies. Our 
main result is that the degeneracy of the first excited state per ground state and per spin is finite 
in the thermodynamic limit. This is very different from the same model on a square lattice where 
a divergence proportional to the linear lattice size is expected. The energy gap for the brickwork 
lattice is obviously 2 J on finite systems and predicted to be the same in the thermodynamic limit. 
Our results suggest that a 2 J gap is universal for planar bimodal Ising spin glasses. The distribution 
of the second contribution to the internal energy has a mode close to zero and we predict that the 
low-temperature specific heat is dominated by the leading term proportional to T~ 2 exp(— 2J/kT). 

PACS numbers: 75.10.Hk, 75.10.Nr, 75.40.Mg, 75.60.Ch 



I. INTRODUCTION 

Although bimodal planar Ising spin glass models are 
very simple in concept, they are extremely controversial. 
One main reason why concerns the energy gap between 
the ground and first excited states. To date, most work 
has been reported for the square lattice where there is 
wide acceptance for the scenario of a critical point only at 
zero temperature 1,2 . In the absence of any contradictory 
evidence or suggestion, we assume this to be the case for 
other planar models, including the brickwork lattice. 

Bimodal models have bond (nearest-neighbor) interac- 
tions of fixed magnitude J and random sign. Both the 
ground and excited states have a very large degeneracy. 
For an infinite square lattice, without open boundaries, 
it is clear that any finite number of spin flips must either 
result in another ground state or an excited state with 
an increase in energy not less than 4 J. For the brickwork 
lattice the corresponding energy gap is clearly 2 J, since 
the lattice coordination number is three. 

About 20 years ago, Wang and Swendsen 3 published 
evidence that the energy gap for the square lattice in the 
thermodynamic limit was 2J. This flew in the face of 
the naive expectation of 4J. Essentially, the claim was 
that it is possible for an infinite number of spin flips to 
provide an excited state with energy only 2J above a 
ground state. The issue here is the noncommutativity 
of the zero-temperature and thermodynamic limits. The 
thermodynamic limit has to be taken first. Nevertheless, 
for the brickwork lattice we do not expect this to be 
an issue and the main interest of this paper is to show 
evidence that this is the case. Both models have the 
same energy gap in the thermodynamic limit, although 
the reasons why are quite different. 

For the square lattice there are three scenarios in 
the literature regarding the energy gap. First, sup- 
port for the 2 J energy gap includes work at finite 
temperatures involving exact computations of parti- 
tion functions^, a worm algorithm^ and Monte Carlo 
simulation 6 . Distributions of excited-state degeneracies 



at arbitrary temperature^ also indicate a 2J gap. Es- 
sentially, it was shown that the degeneracy of the first 
excited state per ground state and per spin diverges in 
the thermodynamic limit. In consequence the 4J gap of 
the finite system is reduced to 2 J. 

Second, Saul and Kardar— reported that the energy 
gap should be 4J as suggested by simple analysis. The 
third published scenario 9 basically claims that the energy 
gap approaches zero in the thermodynamic limit leading 
to power law behaviour for the specific heat. This pos- 
sibility has been discussed at some length in Ref. [Tol . 
although clear conclusions remain unavailable due to dif- 
ficulties related to finite lattice size and extrapolation to 
very low temperature. 

The brickwork lattice (Fig. [lja)) studied in this paper 
is logically equivalent to the hexagonal, or honeycomb, 
lattice (Fig. [U^b)). Very little work has been published to 
date, especially concerning the ground state. The ground 
state energy per bond has been quoted-^ as —0.82 J. The 
entropy per spin has been reported by Aromsawa^ as 
0.02827(5)fc with an energy— 1 .2403 (2) J per spin, in 
good agreement with Ref. [ill The spin glass phase is 
thought to exist for a concentration of negative bondsi^ 
above about 0.067. Work at finite temperatur o 14 ' 15 ! 16 
places the multicritical, or Nishimori, point at the same 
concentration; in agreement if reentrant phase bound- 
aries are absent. 




FIG. 1: (a) The brickwork lattice, (b) The equivalent hexag- 
onal lattice. 

Our calculations of the degeneracies of excited states 
for the brickwork lattice are exact. The temperature 



2 



is fixed and arbitrarily low; we do not use any numer- 
ical value. The lattice is constructed by taking a square 
lattice and diluting bonds in a regular manner to leave 
plaquettes with six perimeter bonds; logically equivalent 
to hexagons as shown in Fig. [TJ The disorder realiza- 
tions are independently quenched random configurations 
of negative bonds in a patch that contains all the frus- 
trated plaquettes. Periodic boundary conditions are used 
in one dimension. The number of sites L for this dimen- 
sion is necessarily even. The cylindrically wound frus- 
trated patch is embedded in an infinite unfrustrated en- 
vironment in the second dimension. 



II. FORMALISM 

An algorithm based on the Pfaffian method^ and de- 
generate state perturbation theor y 18 ! 19 ' 20 for the square 
lattice was adapted to evaluate the degeneracies of ex- 
cited states for the brickwork lattice. The main points 
of this procedure are given in the following. From the 
square lattice, we dilute bonds in a regular manner to 
define the brickwork lattice. Using the fermion decora- 
tion of bonds (one either side) , a brickwork plaquette has 
eight fermions inside (filled circles) and six others across 
the bonds as shown in Fig. [2j 
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(a) (b) 

FIG. 2: (a) A square plaquette as in Ref. [20. (b) A brickwork 
plaquette obtained by dilution of the central bond (denoted 
by a dotted line) between two square plaquettes. 

As for the bimodal Ising spin glass on the square lat- 
tice, the partition function is given by Z ~ (det-D) 1 / 2 
where D is a, real skew-symmetric AN x AN matrix for a 
lattice with TV sites. D is a singular matrix at zero tem- 
perature with zero eigenvalues which are equal in number 
to the number of frustrated plaquettes. Each eigenvalue 
approaches zero according to the form 

e = ±\x exp(-2 Jr/kT), (1) 

where r is an integer (an order of perturbation theory) 
and A is a real number. The quantities r and X can be 
exactly evaluated by degenerate state pertubation the- 
ory. The ground-state energy and entropy of the system 
can be defined similarly as for the square lattice^. It is 
equivalent to expressing the ground-state degeneracy as 

M = Y[X. (2) 



where the product is over all the positive defect eigenval- 
ues. 

The gauge-invariant method is applied similarly as for 
the square lattice^ 9 - to separate the singularities of D 
for the brickwork lattice using real matrices as follows. 
At zero temperature, the perfect system (no frustra- 
tion) Green's functioni 9 - go is obtained by transforming 
D into a plane-wave basis and inverting an 8 x 8 ma- 
trix. The nonzero elements of go are only across bonds 
and localized inside plaquettes. Across bonds we have 
.9o (+i — ) = — 9o{~ j +) = g w here the matrix indices are 
defined in Fig. 4 of Ref. T9j. Within a plaquette go is as 
given by the following matrix^: 
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(3) 



where the elements of go are with respect to the bond 
basis of a plaquette as shown in Fig. [3] 
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FIG. 3: The states in the bond basis of a brickwork plaquette. 
It is convenient to use the bond basis introduced in Ref. 0, 
that is | ±) = ;j =(| a ) ± C | P)) where | a) and | /3) are shown 
in Fig. 4 of Ref. and £ represents the sign of the bond. 

We set the bimodal problem for the brickwork lattice 
by placing negative bonds at random. To reduce the 
complexity, we perform gauge transformations as in Ref. 
Il9l so that the negative (defect) bonds are vertical only. 
We classify plaquettes into four types as shown in Fig. 2] 



(a) (b) 



(c) (d) 

FIG. 4: Heavy lines denote negative (defect) bonds. There are 
four types of possible plaquettes after gauge transformation to 
vertical defect bonds, (a) and (b) are unfrustrated plaquettes 
while (c) and (d) are frustrated plaquettes. 
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To determine ground state properties, degenerate state 
perturbation theory is applied at arbitrarily small finite 
temperatures. We write exactly D = Dq + 5D\ where 
(5 = 1 — t, with t = tanh(J/feT), is used as a parameter 
for a perturbation expansion. The matrix Dq is singular 
and has defect eigenvectors | d) corresponding to defect 
(zero) eigenvalue, that is Dq \ d) = 0, localized on each 
frustrated plaquette, similar to those found for the square 
lattice^. 

At first order the matrix D\ , which is 2 x 2 block diag- 
onal across bonds, is diagonalized in the basis of the de- 
fect eigenvectors. We also require the continuum Green's 
function G c = (1 - P)g c {l - P) where P = £ J d)(d\ 
and g c = go + goUgo with U defined similarly to Ref. ll9L 
We can write g c = g c \ +g C 2 where g c \ has matrix elements 
connecting the basis states within a brickwork plaquette 
while g C 2 connects only basis states across a single bond, 
that is g C 2{+,— ) = — <?c2(— ,+) = \ and is only rele- 
vant for excited states. Although g can take us to the 
fcrmions associated with diluted bonds (labeled 7 — 8 in 
Fig. [3]), it can be proven that these matrix elements do 
not effect the calculation of the partition function since 
Di across that bond is equal to zero (there is no energy) ; 
we can disregard them. An alternative Pffaffian based 
on three nodes per site has been given in Ref. [13 but it 
cannot be adapted to our defect problem. 

The matrix g c i can be considered in the subbasis with 
only six fcrmions (labeled 1 — 6 in Fig. [3J without any 
change of the gauge-invariant ground or excited state 
properties. We can also arrange to have g c \ orthogo- 
nal to defect states, that is g c \ \ d) = 0, by under- 
standing that we can add any matrix A to g c \ as long 
as (1 — P)A(1 — P) =0. This reduces the number of 
arithmetic operations for the calculation of excitations. 
The matrix g c \ can be presented for an unfrustrated pla- 
quette as 



9cl 



-1 1 -1 s s 

1 1 -1 s s 
-1 -1 -1 s s 

1 1 1 s s 

— s — s — s —s 1 

— s — s — s —s —1 



(4) 



where s = 1 for an unfrustrated plaquette with no neg- 
ative bond (Fig. HJa)) and s = — 1 for an unfrustrated 
plaquette with two negative bonds (Fig. BJb)). The ma- 
trix U only occurs for plaquettes with two defect bonds. 
In detail, C/34 = — E/43 = —2. The matrix g c i for a frus- 
trated plaquette is given by 



9d 



- 


-2 


2 


-1 


— s 





2 





1 


-2 





s 


-2 


-1 








-2s 


— s 


1 


2 








s 


2s 


s 





2s 


— s 





2 





— s 


s 


-2s 


-2 






(5) 



where s = 1 for a frustrated plaquette with a left negative 
bond (Fig. He)) an d s = — 1 otherwise (Fig. H£d)). The 



corresponding defect states are 
1 



d) 



V6 



(| 1)+ I 2)+ I 3)+ I 4) -s (| 5)+ I 6))) (6) 



The prefactor in Eq. (p} is essentially determined by the 
normalization of this vector. 

At second order, we diagonalize D2 = D\g c \D\. To 
continue for higher orders, we require the Green's func- 
tions G ri as given in Ref. [7J, obtained from previous 
orders; that is for states whose degeneracy has already 
been lifted. The general rule for D r (at rth order) 
can be expressed as D r = D r _i(l + G r _2-Dr-2) • • ■ (1 + 

The calculation of the internal energy and the specific 
heat for the brickwork lattice is similar to the square 
lattice as described in Ref. 0- The internal energy is 
given by 



U=J2 e- 2Jm ' kT U n 



(7) 



where Uq is the ground state energy and U m — 
-2 m JTri? m , for m > 0, and 



R = D l9cl (1 + Did) (1 + D 2 G 2 ) . . . (1 + D r 



<G rma J. 
(8) 

r max is the highest order of pertubation theory required. 

However, there is one essential distinction between the 
square and brickwork lattices. Since we need three colors 
to color the brickwork lattice, this means that the color 
rules described in Ref. are invalid and it follows that 
U m 7^ for all m. The specific heat per spin can be 
expressed in terms of the internal energy as 



1 dU _ 1 / 2J \ ^ 

v / in— 1 



-2Jm/kT T r 



(9) 



The degeneracy of the ith excited state is given as M, 

••Mq 



l/Mx \2\ 



III. RESULTS 

Fig. E] shows the distributions of jj^* for the spin glass 
with (£+1) x L spins. It is clear that the most likely value 
scales as the number of spins. We also consider a different 
shape of the patch boundary by changing the lattice size 
to (2L + 1) x L so that the equivalent hexagonal lattice 
has a more balanced shape. The distribution of ^ with 
(2L + 1) x L spins also shows again a most likely value 
that scales as the number of spins as shown in Fig. [6] 
Also, the value is roughly the same. In both cases, it is 
clear that the leading term of the specific heat grows like 
c v ~ T~ 2 exp(-2 J/kT) indicating a 2J excitation gap. 
Moreover, the distributions of are extreme, do not 
self-average and are neither normal or lognormal. This 
is consistent with the square lattice as in Ref. 0- 
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FIG. 5: (Color online) Distributions for |p. 



0.8 



with var- 



M (L+l)xL 

ious values of lattice size L with (L + 1) xi spins. Each dis- 
tribution includes 10 5 disorder realizations, except for L = 96 
which has 30, 000. 
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FIG. 6: (Color online) Distributions for lf , 
various values of lattice size L with (2L + 1) x L spins. Each 
distribution includes 10 5 disorder realizations. 



The distributions of j^- clearly show extreme-value 
distributed properties with a fat tail. It is thus reason- 
able to analyze the distributions according to the Frechet 
distribution with the probability density function: 



/' 







exp 







( 10 ) 

where the parameters /i, p and £ indicate location, shape 
and scale of the distribution respectively. The mode can 



be calculated by x — [i + (3 



We estimate the 



parameters by a maximum likelihood estimator 21 to fit 
actual disorder realizations. It is found that the Frechet 
distribution cannot fit exactly the peak of our actual dis- 
tribution; the quality of the fit is quite poor and also 
deteriorates with respect to increasing L. An example 
for L = 48 with (2L + 1) x L spins is shown in Fig. [TJa). 



We have also tried to polish the fit using the Levenberg- 
Marquardt method 22 to fit the bin data, setting initial 
values of parameters equal to the values obtained from 
the algorithm of Hosking 21 . This provides alternative pa- 
rameters and the curve is shown in Fig. [7](b) . Although 
the quality of the fit looks a little better and the position 
of the mode somewhat improves, we are not convinced 
that Eq. (TIT)]) should necessarily describe the distribu- 
tion exactly. For comparison, at different values of L, 



FIG. 7: Distributions for 



with L 



48 and 



M (2L+l)xL 

(2L + 1) X L spins, (a) The line uses the parameters obtained 
from the algorithm of Hosking— fitted to actual disorder real- 
izations, (b) The line uses the parameters obtained from the 
algorithm of Levenberg-Marquardt— fitted to bin data using 
the Hosking parameters as initial values. 



the fitted distributions divided by f(x) are shown in Fig. 
isJU It is also useful to present the mode of the fitted dis- 
tributions as a function of L in Fig. [9] We trust that the 
/mode converges; there is no reason to believe otherwise. 



We have also found distributions for the second con- 
tribution to the internal energy, ^ — \{^p) 2 - These 
are shown in Figs. [10] and [11] The most fikely value 
is very close to zero. The skewness also suggests the 
dominance of the first excitations. We thus believe that 
higher excitations are unlikely to change our conclusion 
that the energy gap is 2 J. The mode of ^ alone scales 
as ((L + 1)L) 2 and ((2L+1)L) 2 as shown in Figs. [T2land 
[T51 respectively. 
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FIG. 8: (Color online) Frechet distributions from fitting the 
bin data in Fig. [6] by the algorithm of Levenberg-Marquardt 
using the Hosking parameters as initial values. 
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FIG. 10: (Color online) Distributions for the second term in 
the specific heat with various values of lattice size L with 
(L + 1) X L spins. Each includes 10 s disorder realizations, 
except for L = 96 which has 30, 000. 
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FIG. 9: The mode of the fitted distributions with various 
values of lattice size L with (2L + 1) x L spins. 
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FIG. 11: (Color online) Distributions for the second term in 
the specific heat with various values of lattice size L with 
(2L + 1) X L spins. Each includes 10 5 disorder realizations. 



IV. CONCLUSIONS 

In conclusion, we have reported exact results for the 
excitations of the bimodal Ising spin glass on the brick- 
work lattice by expanding in arbitrary temperature from 
the ground state. This is complimentary to the more 
usual extrapolation from finite temperature. 

We find that the energy gap is 2 J for both finite and in- 
finite lattices. The thermodynamic limit is trivial in con- 
trast to the difficulties associated with the square lattice. 
Our result may suggest that a 2 J energy gap is universal 
for planar bimodal Ising spin glasses in the thermody- 
namic limit. For instance, the triangular lattice could 
very well be expected to behave similarly to the square 
lattice since its plaquettes can be colored using just two 
colors; the brickwork lattice requires three colors. 

As a final note, we expect a correlation length £ ~ 
exp(2J//cT) in probable agreement^ ^ 23 ' 24 with the 
square lattice. Our reasoning is based on the construction 



of correlation functions using reciprocal defect a ' 18 > 20 
and closed polygons. The essential point is that, for a 
finite lattice, the correlation functions must be analytical 
functions of t and thus of S ~ exp(— 2J/fcT). Compar- 
ison with the asymptotic expression exp(—R/£)/RP for 
correlation functions at large separation i?, allows us to 
deduce that £ _1 is also an analytical function of 6. If the 
thermodynamic limit is trivial for the energy gap, then it 
is very likely to be so for the correlation length. It is also 
known that the fully frustrated brickwork lattice has a 
constant correlation length 25 and that the ground state 
is not critical, unlike the Villain model 2 ^ that has a non- 
analytical free energy although ~ 6 nevertheless 2 ^. 
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FIG. 12: (Color online) Distributions for ffi( (L+ \ )xL ) 2 with 
various values of lattice size L with (L + 1) x L spins. Each 
distribution includes 10 5 disorder realizations, except for L — 
96 which has 30, 000. 
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FIG. 13: (Color online) Distributions for ^1 
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various values of lattice size L with (2L + 1) x L spins. Each 
distribution includes 10 5 disorder realizations. 
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